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Abstract 

If supersymmetry is discovered at the LHC, the measured spectrum of superpartner 
masses and couplings will allow us to probe the origins of supersymmetry breaking. 
However, to connect the collider-scale Lagrangian soft parameters to the more funda- 
mental theory from which they arise, it is usually necessary to evolve them to higher 
scales. The apparent unification of gauge couplings restricts the possible forms of new 
physics above the electroweak scale, and suggests that such an extrapolation is pos- 
sible. Even so, this task is complicated if the low scale spectrum is not measured 
completely or precisely, or if there is new physics at heavy scales beyond the reach of 
collider experiments. In this work we study some of these obstacles to running up, and 
we investigate how to get around them. Our main conclusion is that even though such 
obstacles can make it very difficult to accurately determine the values of all the soft 
parameters at the high scale, there exist a number of special combinations of the soft 
parameters that can avoid these difficulties. We also present a systematic application 
of our techniques in an explicit example. 
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1 Introduction 



It is a remarkable feature of the minimal supersymmetric extension of the standard model (MSSM) 
that the SU(3) C , SU(2)l, and U(l)y gauge couplings unify at a very high scale, of order 
10 16 GeV [UI2]- Furthermore, the matter fields of the MSSM, with the exception of the pair 
of Higgs doublets, have precisely the quantum numbers of three sets of 5 © 10 represen- 
tations of SU(5). These properties may be purely accidental, but they do suggest a more 
symmetric unified structure at energies only slightly below the Planck scale. They also offer 
a tantalizing hint of the structure of Nature at scales well above what we will be able to 
probe directly with colliders such as the LHC. 

Taking these features as being more than accidental, we obtain significant constraints 
on the types of new physics that can arise between the electroweak and the grand unifi- 
cation (GUT) scales. Any new phenomenon that enters the effective theory in this energy 
range ought to maintain the unification of couplings, and should be consistent with a (pos- 
sibly generalized) GUT interpretation. The simplest scenario is a grand desert, in which 
there is essentially no new physics at all below the unification scale Mqut- In this case, 
if supersymmetry is discovered at the Tevatron or the LHC, it will be possible to extrap- 
olate the measured soft supersymmetry breaking parameters to much higher scales using 
the renormalization group (RG). Doing so may help to reveal the details of supersymmetry 
breaking, and possibly also the fundamental theory underlying it. 

If supersymmetry is observed in a collider experiment, it will be challenging to extract 
all the supersymmetry breaking parameters from the collider signals. While some work has 
been put into solving this problem [3], there is still a great deal more that needs to be 
done. The parameters extracted in this way will be subject to experimental uncertainties, 
especially if the supersymmetric spectrum is relatively heavy. There will also be theoretical 
uncertainties from higher loop corrections in relating the physical masses to their running 
values [3]. These uncertainties in the supersymmetry breaking parameters, as well as those 
in the supersymmetric parameters, will complicate the extrapolation of the soft masses to 
high energies [5J E|. Much of the previous work along these lines has focused on running the 
parameters of particular models from the high scale down. This is useful only if the new 
physics found resembles one of the examples studied. Our goal is to study the running from 
low to high [TJ. 

Evolving the soft parameters from collider energies up to much higher scales can also be 
complicated by new physics at intermediate energies below Mqut- The apparent unification 
of gauge couplings suggests that if this new physics is charged under the MSSM gauge 
group, it should come in the form of complete GUT multiplets or gauge singlets. Indeed, the 
observation of very small neutrino masses already suggests the existence of new physics in 
the form of very heavy gauge singlet neutrinos [8]. With new physics that is much heavier 
than the electroweak scale, it is often very difficult to study it experimentally, or to even 
deduce its existence. If we extrapolate the MSSM parameters without including the effects 
of heavy new physics, we will obtain misleading and incorrect values for the high scale values 
of these parameters [9]. 
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In the present work we study some of these potential obstacles to the RG evolution of 
the MSSM soft parameters. In Section [2] we investigate how uncertainties in the low-scale 
parameter values can drastically modify the extrapolated high-scale values. We focus on 
the so-called S term (a.k.a. the hypercharge D term) within the MSSM, which depends on 
all the soft masses in the theory, and can have a particularly large effect on the running if 
some of these soft masses go unmeasured at the LHC. In Sections [3] and HI we study two 
possible examples of heavy new physics. Section [3] investigates the effects of adding complete 
vector-like GUT multiplets on the running of the soft parameters. Section H] describes how 
including heavy Majorana neutrinos to generate small neutrino masses can alter the running 
of the MSSM soft parameters. In Section [5] we combine our findings and illustrate how they 
may be put to use with an explicit example. Finally, Section[6]is reserved for our conclusions. 
A summary of some useful combinations of scalar soft masses is given in an Appendix. 

Our main result is that the high scale values of many Lagrangian parameters can be very 
sensitive to uncertainties in their low-scale values, or to the presence of heavy new physics. 
However, in the cases studied we also find that there are particular combinations of the 
Lagrangian parameters that are stable under the RG evolution, or that are unaffected by 
the new physics. These special parameter combinations are therefore especially useful for 
making a comparison with possible high-scale theories. 

Throughout our analysis, we simplify the RG equations by setting all flavor non-diagonal 
soft terms to zero and keeping only the (diagonal) Yukawa couplings of the third genera- 
tion. Under this approximation, we work to two loop order for the running of the MSSM 
parameters, and interface with Suspect 2.3.4 [10] to compute one-loop threshold correc- 
tions at the low scale. For concreteness, we take this scale to be 500 GeV. The additional 
running due to new physics introduced at scales much larger than the weak scale is only 
performed at one-loop. We also implicitly assume that the mass scale of the messengers that 
communicate supersymmetry breaking to the visible sector lies at or above the GUT scale, 
Mqut — 2.5 x 10 16 GeV. Even so, our methods and general analysis will also be applicable 
to scenarios that have lighter messenger particles, such as gauge mediation [TH [12] . We also 
neglect the effects of hidden sector running, which can be significant if there are interacting 
states in the hidden sector significantly lighter than Mqut [EH [EH [32] . While this manuscript 
was in preparation, methods similar to those considered in the present work were proposed 
in Ref. [15] to deal with these additional uncertainties in the high scale values of the soft 
parameters. 

In this work, we focus on low-energy supersymmetric models, and particular forms of in- 
termediate scale new physics. Despite this restriction, we expect that our general techniques 
will be applicable to other solutions of the gauge hierarchy problem, or to more exotic forms 
of new intermediate scale physics. 
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2 Uncertainties Due to the S Term 



The one-loop renormalization group (RG) equations of the MSSM soft scalar masses have 
the form [16] 

(16^)^=^- Z9 2 a CnM a f + 6 -gtY t S, (1) 

*' a=l,2,3 ° 

where t = \n(Q/Mz), Xi is a function of the soft squared masses and the trilinear couplings, 
M a denotes the a-th gaugino mass, and the 5* term is given by 

S = Tr(Ym 2 ) (2) 

2 2ij./2 o2, 2, 2 2\ 

= m H u - m H d + tr (m Q -2m u + m E + m D - m L ) 

where the first trace runs over all hypercharge representations, and the second runs only 
over flavors. 

The S term is unusual in that it connects the running of any single soft mass to the soft 
masses of every other field with non-zero hypercharge. Taking linear combinations of the 
RG equations for the soft scalar masses, the one-loop running of S in the MSSM is given by 

(16vr 2 )^ = -2 6 1 ^5, (3) 

where bi = —33/5 is the one- loop beta-function coefficient. Using Eq. ([3]) in Eq. ([1]) and 
neglecting the Yukawa-dependent terms Xi (which are expected to be small for the first and 
second generations) the effect of S ^ is to shift the high scale value the soft mass would 
have had were S(to) = by an amount 



Tr(Y* 



gfCQ 
Ui(*o) 



- 1 



S(t ). (4) 



The one-loop RG equation for S is homogeneous. Thus, if S vanishes at any one scale, it 
will vanish at all scales (at one-loop). In both minimal supergravity (mSUGRA) and simple 
gauge-mediated models, 5* does indeed vanish at the (high) input scale, and for this reason 
the effects of this term are often ignored. 

From the low-energy perspective, there is no reason for S(to) to vanish, and in many 
cases its effects can be extremely important. Since g\ grows with increasing energy, the 
mass shift due to the S term grows as well. For t ~ tu z and t — > tour, the prefactor in 
Eq. (BJ is about (0.13) Fj. The value of S(t ) depends on all the scalar soft masses, and 
the experimental uncertainty in its value will therefore be set by the least well-measured 
scalar mass. In particular, if one or more of the soft masses aren't measured at all, S(to) 
is unbounded other than by considerations of naturalness. Fortunately, this uncertainty 
only affects the soft scalar masses. The S term does not enter directly into the running of 
the other soft parameters until three- loop order [HI [17] , and therefore its effects on these 
parameters is expected to be mild. 
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There is also a theoretical uncertainty induced by the S term. Such a term is effectively 
equivalent to a Fayet-Iliopoulos (FI) D term for hypercharge |17j . To see how this comes 
about, consider the hypercharge D-term potential including a FI term, 



miWW (5) 

where we have also included the soft scalar masses. Eliminating the Di through its equation 
of motion, we find 

2 = ■ ■ ■ - \e - f (y,4>iQiA - x> + \j\^ Y ^ ) ^ ( 6 ) 

Thus, except for the constant addition to the vacuum energy, the effect of the FI term is to 
shift each of the soft squared masses by an amount 

m \ ■= ™\ + y| gi Y i f • (7) 

The low-energy observable quantities are the {mj}, not the {fnf\- Since we can't extract 
the shift in the vacuum energy, the low-energy effects of the FI term are therefore invisible 
to us, other than the shift in the soft masses. This shift is exactly the same as the shift due 
to the S term. Let us also mention that the S term, as we have defined it in Eq. (j2J), runs 
inhomogeneously at two-loops and above, so the exact correspondence between a hypercharge 
D term and the S term of Eq. (j2J) does not hold beyond one-loop order. 

A simple way to avoid both the large RG uncertainties in the soft masses and the theoret- 
ical ambiguity due to the S term is apparent from Eq. PJ . Instead of looking at individual 
soft masses, it is safer to consider the mass differences 

Yjml-Yiml (8) 

for any pair of fields. These differences are not affected by the mass shifts of Eq. They 
are also independent of the value of the FI term. 

In the rest of this section, we show how the S term can complicate the running of the 
soft masses to high energies with a particular example. If one of the scalar soft masses is 
unmeasured, it is essential to use the linear mass combinations given in Eq. (jSJ) instead of 
the individual masses themselves. We also discuss how the special RG properties of the S 
term provide a useful probe of high scale physics if all the scalar soft masses are determined 
experimentally. 

2.1 Example: SPS-5 with an Unmeasured Higgs Soft Mass 

To illustrate the potential high-scale uncertainties in the RG-evolved soft parameters due to 
the S term, we study the sample mSUGRA point SPS-5 [18] under the assumption that one 
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Figure 1: Deviations in the running of some of the SPS-5 soft masses due to setting m 2 Hd = 
(1000 GeV) 2 at the low scale. The solid lines show the running of mn 2 ^, fn 2 E , and m\ with this 
perturbation, while the dashed lines show the unperturbed running of these soft masses. The 
unperturbed low-scale value of the down-Higgs soft mass is fn 2 Hd ~ (235 GeV) 2 . 

of the scalar soft masses goes unmeasured at the LHC. If this is the case, the S parameter 
is undetermined, and the precise values of the high-scale soft terms are no longer precisely 
calculable. Even if the value of the S term is bounded by considerations of naturalness, the 
uncertainties in the high-scale values of the soft scalar masses can be significant. 

The SPS-5 point is defined by the mSUGRA input values m = 150 GeV, rayi = 
300 GeV, A = —1000, tan/3 = 5, and sgn(fj,) > 0, at Mqut- The mass spectrum for 
this point has relatively light sleptons around 200 GeV, and somewhat heavier squarks with 
masses near 400-600 GeV. The LSP of the model is a mostly Bino neutralino, with mass 
close to 120 GeV. The perturbation we consider for this point is a shift in the down-Higgs 
soft mass, m 2 H . 

Of the soft supersymmetry breaking parameters in the MSSM, the soft terms associated 
with the Higgs sector can be particularly difficult to deduce from LHC measurements. At 
tree-level, the independent Lagrangian parameters relevant to this sector are pQ 

v, tan/5, [A, Ma, (9) 

where v ~ 174 GeV is the electroweak breaking scale, tan/5 = v u /vd is the ratio of the H u 
and H d VEVs, \i is the supersymmetric //-term, and Ma is the pseudoscalar Higgs boson 
mass. Other Higgs-sector Lagrangian parameters, such as m 2 H and m 2 H , can be expressed 
in terms of these using the conditions for electroweak symmetry breaking in the MSSM. 
Among the Higgs sector parameters listed in Eq. ([9]), only the value of v is known. The 
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Figure 2: Deviations in the running of some of the SPS-5 soft masses due to setting m 2 Hd = 
(1000 GeV) 2 at the low scale. The solid lines show the running of 2 , m 2 j i2 , and fn 2 ^ 

with this perturbation, while the dashed lines show the unperturbed running of these soft masses. 

value of /i can potentially be studied independently of the Higgs scalar sector by measuring 
neutralino and chargino masses and couplings 0, HH EH] , although it is likely to be poorly 
determined if only hadron colliders are available. A number of observables outside the Higgs 
sector may also be sensitive to tan /3, especially if it is large, tan (3 > 20. For example, 
the dilepton invariant mass distributions in the inclusive 2£ + jets + IpT channel can vary 
significantly depending on the value of tan (3, but this dependence is such that the value of 
tan (3 can at best only be confined to within a fairly broad ranges [2U] . Determining Ma at 
the LHC typically requires the discovery of one of the heavier MSSM Higgs boson states. 
Finding these states can also help to determine tan/? [21]. Unfortunately, the LHC reach 
for the heavier Higgs states is limited, especially for larger values of Ma and intermediate or 
smaller values of tan/3 < 20 [22]. If none of the heavier Higgs bosons are found at the LHC, 
it does not appear to be possible to determine both of the Higgs soft masses, rn 2 Hu and rn 2 Hd . 

For the low-scale parameters derived from SPS-5, the pseudoscalar Higgs mass Ma is 
about 700 GeV. With such a large value of Ma, and tan/3 = 5, only the lightest SM-like 
Higgs boson is within the reach of the LHC [22]. Motivated by this observation, we examine 
the effect of changing the low-scale value of m 2 ^ on the running of the other soft parameters. 
The actual low-scale value of m 2 ^ is about (235 GeV) 2 . The perturbation we consider is to 
set this value to (1000 GeV) 2 , while keeping tan/3 fixed. Such a perturbation does not ruin 
electroweak symmetry breaking, and tends to push the heavier higgs masses to even larger 
values. In the present case, the heavier Higgs masses increase from about 700 GeV to over 
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1200 GeV. 

The effects of this perturbation in rn? Hd on some of the soft scalar masses are shown 
in Figs. [H and El In these plots, we show rrii = mf / ^J\mf\. The deviations in the soft 
masses are substantial, and the S parameter is the source of this uncertainty. In addition 
to the S term, varying m 2 Hd can also modify the running of the Higgs mass parameters and 
the down-type squarks and sleptons through the Yukawa- dependent terms X, in the RG 
equations, Eq. (CQ). In the present case these Yukawa- dependent effects are very mild since 
for tan (3 — 5, the b and r Yukawas are still quite small. This can be seen by noting the 
small difference between the perturbed running of m 2 Jl 2 and mfj 3 in Fig. [2j For larger values 
of tan /3, the down-type Yukawa couplings can be enhanced and this non-S" effect from rn? Hd 
can be significant. However, we also note that as these Yukawa couplings grow larger, it 
becomes much easier for the LHC to detect one or more of the heavier Higgs states. To the 
extent that the Yukawa-dependent shifts can be neglected, the linear mass combinations of 
Eq. (jSj) remove most of the uncertainty due to an unknown m 2 Hd in the running of the soft 
masses that are measured. The effect of not knowing rn 2 Hd has only a very small effect on 
the running of the gaugino masses and the trilinear terms. 

In this example we have assumed that m? H is the only unmeasured soft scalar mass. 
Several of the other soft scalar masses may be difficult to reconstruct from LHC data as 
well. For example, within many SUSY scenarios the third generation squarks and some of 
the heavier sleptons have very small LHC production rates. If there are other unmeasured 
soft scalar masses besides m? H , the uncertainties due to the S term in the extrapolation 
of the measured scalar soft masses will be even greater than what we have presented here. 
The mass combinations of Eq. (jSJ) will be necessary to study the high scale supersymmetry 
breaking spectrum in this case. 



2.2 Origins and Uses of the S Term 

While the S term can complicate the extrapolation of the soft scalar masses if one of them 
goes unmeasured, the simple scale dependence of this term also makes it a useful probe of the 
high scale theory if all the masses are determined. The essential feature is the homogeneous 
RG evolution of the S term, given in Eq. ([3]), which is related to the non-renormalization of 
FI terms in the absence of supersymmetry breaking. 

A non- vanishing S term can arise from a genuine FI term present in the high scale theory, 
or from non-universal scalar soft masses at the high input scale. The size of a fundamental 
hypercharge FI term, £, is naively on the order of the large input scale. Such a large value 
would either destabilize the gauge hierarchy or lead to U(l)y breaking at high energies. 
However, the non-renormalization of FI terms implies that it is technically natural for £ 
to take on much smaller values. In this regard, a small value for £ is analogous to the /x 
problem in the MSSM. Adding such a FI term to mSUGRA provides a simple one-parameter 
extension of this model, and can have interesting effects [23j 121] . Note, however, that in a 

1 The values of fi and Bji also change, although the variation in /j, is very mild: /i ~ 640 GeV — > 670 GeV. 



GUT where U(l)y is embedded into a simple group, a fundamental hypercharge FI term in 
the full theory is forbidden by gauge invariance. 

It is perhaps more natural to have the S term emerge from non-universal scalar soft 
masses j25j EE]. This is true even in a GUT where U(l)y is embedded into a simple group. 
Within such GUTs, the contribution to S from complete GUT multiplets vanishes. However, 
non-zero contributions to S can arise from multiplets that are split in the process of GUT 
breaking. For example, in 577(5) with H u and H d states embedded in 5 and 5 multiplets, 
a non-zero low-energy value of S can be generated when the heavy triplet states decouple 
provided the soft masses of the respective multiplets are unequal. Whether it is zero or 
not, the low-scale value of the S term provides a useful constraint on the details of a GUT 
interpretation of the theory. 

So far we have only considered the S term corresponding to U(l)y. If there are other 
gauged U(l) symmetries, there will be additional S term-like factors for these too. In fact, 
the homogeneity of the S term evolution also has a useful implication for any non-anomalous 
global £7(1) symmetry in the theory. The only candidate in the MSSM is U(1)b-l, up to 
linear combinations with U(l)y [27]. Let us define Sb-l by the combination 

Sb-l = Tr(Q B - L m 2 ) (10) 
= tr{2m 2 Q — vrt^j — m 2 D — 2m 2 L + m 2 E ), 

where the second trace runs only over flavors|§ We can think of Sb-l as the effective S term 
for a gauged U(1)b-l m the limit of vanishing coupling. At one loop order, the RG running 
of Sb-l is given by 

{167[2) ^r = l Tr (QB-LY)giS (11) 
= n g — gi S, 

where n g = 3 is the number of generations. If S is measured and vanishes, Sb-l provides 
a second useful combination of masses that is invariant under RG evolutionjl and yields an 
additional constraint on possible GUT embeddings of the theory. 



3 New Physics: Complete GUT Multiplets 

As a second line of investigation, we consider the effects of some possible types of new 
intermediate scale physics on the running of the MSSM soft terms. If this new physics is 
associated with supersymmetry breaking as in gauge mediation [TT], then it is of particular 
interest in its own right. Indeed, in this case the low-energy spectrum of soft terms may 
point towards the identity of the new physics after RG evolution. On the other hand, there 

2 Up to flavor mixing, we can also define an independent Sb-l f° r each generation. 

3 This non-evolution of mass combinations corresponding to non-anomalous global symmetries persists 
at strong coupling. In models of conformal sequestering, this can be problematic 13J. 
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are many kinds of possible new intermediate scale physics that are not directly related to 
supersymmetry breaking. The existence of this type of new phenomena can make it much 
more difficult to deduce the details of supersymmetry breaking from the low-energy soft 
terms. 

A useful constraint on new physics is gauge coupling unification. To preserve unification, 
the SM gauge charges of the new physics should typically be such that all three MSSM gauge 
beta functions are modified in the same way@ This is automatic if the new matter fills out 
complete multiplets of a simple GUT group into which the MSSM can be embedded. Such 
multiplets can emerge as remnants of GUT symmetry breaking. 

As an example of this sort of new physics, we consider vector-like pairs of complete SU{5) 
multiplets. For such multiplets, it is possible to write a down a supersymmetric mass term 
of the form 

WDjlXX, (12) 

where X and X denote the chiral superfields of the exotic multiplets. We also assume 
that the exotic multiplets have no significant superpotential (Yukawa) couplings with the 
MSSM fields. Under these assumptions, the exotic SU (5) multiplets can develop large masses 
independently of the details of the MSSM. They will interact with the MSSM fields only 
through their gauge interactions. 

If the supersymmetric mass ft is much larger than the electroweak scale, it will be very 
difficult to deduce the presence of the additional GUT multiplets from low-energy data alone. 
Moreover, an extrapolation of the measured soft masses using the RG equations appropriate 
for the MSSM will lead to incorrect values of the high-scale parameters. In this section, we 
characterize the sizes and patterns of the deviations in the high scale soft spectrum induced 
by additional vector-like GUT multiplets. We also show that even though the new matter 
interferes with the running of the MSSM soft parameters, it is often still possible to obtain 
useful information about the input spectrum, such as the relative sizes of the gaugino masses 
and whether there is inter-generational splitting between the soft scalar masses. 



3.1 Shifted Gauge Running 

The main effect of the exotic GUT multiplets is to shift the running of the SU(3) C , SU(2)l, 
and U (l)y gauge couplings. Recall that in the MSSM, the one-loop running of these couplings 
is determined by 

dg~ 2 _ b a 

with {pi, 62, 63) = (—33/5, —1,3). The presence of a massive GUT multiplet shifts each of 
the b a up by an equal amount above the heavy threshold at t — ti — ln(/2/Mz), 

Ab= -iV 5ffi 5-3iV 1 oei-o + ... (14) 



'For an interesting exception, see Ref. [28] . 
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where iV 5e g is the number of additional 5©5 representations and iV 10el - is the number of 
10 © 10's. The modified one-loop solution to the RG equations is therefore, 



( J " I £T 2 (*o) + - t ) + 0(t - t,) t > tj. {lb} 

It follows that the unification scale is not changed, but the value of the unified gauge coupling 
is increased. Note that the number of new multiplets is bounded from above for a given value 
of p, if gauge unification is to be perturbativeEl 

The change in the gauge running shifts the running of all the soft parameters, but the 
greatest effect is seen in the gaugino masses. At one-loop, these evolve according to 

dM a b a 2 

-df = -^9aM a . (16) 

It follows that M a f ' g\ is RG invariant above and below the heavy threshold. If the threshold is 
also supersymmetric, M a will be continuous across it at tree-level. Since g a is also continuous 
across the threshold at tree-level, the addition of heavy vector-like matter does not modify 
the one-loop scale invariance of the ratio M a j g 2 a . This holds true whether or not the new 
matter preserves gauge unification, but it is most useful when unification holds. When 
it does, the measurement of the low-energy gaugino masses immediately translates into a 
knowledge of their ratio at M GUT [29] . 

From Eq. (IT51) and the one-loop scale invariance of M a /g 2 , the shift in the gaugino masses 
due to the additional matter is 



M„(i) = MJt) 



57T 



1 



(17) 



where g a and M a denote the values these parameters would have for Ab = (i.e. the values 
obtained using the MSSM RG equations). For t = tour, the shift coefficient is identical for 
a = 1,2,3 provided gauge unification occurs. The shift in the running of the gauge couplings 
and the gaugino masses due to seven sets of 5 © 5's at 10 11 GeV is illustrated in Fig. [3J 
An unperturbed universal gaugino mass of Mi/ 2 = 700 GeV is assumed. Both the values of 
the unified gauge coupling and the universal gaugino mass at Mqut are increased by the 
additional multiplets. 

The running of the soft masses also depends on the running of the gauge couplings and 
gaugino masses, and is modified by the appearance of new matter. At one-loop order, in the 
limit that we can neglect the Yukawa couplings, it is not hard to find the shifts in the soft 
masses. We can divide these shifts into two contributions, 

mlit) =m 2 i + Am 2 lx + Am l> ( 18 ) 



5 Note that the parts of \x corresponding to the doublet and triplet components of the 5 will run differently 
between Mqut and the intermediate scale. This will induce additional threshold corrections that we do not 
include in our one-loop analysis. 
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Figure 3: The shift in the running of the gauge couplings (left) and the gaugino masses (right) due 
to 7 sets of 5 © 5's with mass Jl = 10 11 GeV. The universal gaugino mass is taken to be 700 GeV. 



where m? is the value of the soft mass obtained by running the measured value up in the 
absence of the new matter, Amf x is the shift due to the modified gaugino masses, and Amf s 
is due to the change in the running of the S term. 

The first shift, Am^, can be obtained straightforwardly from Eqs. (jTHl3|H6p . and is given 

by 



E 2 ^ 



M„ 



9a 



AL 



(19) 



where Cf the Casimir invariant of the a-th gauge group and 



AL 



JJ; / 4 _ 4 \ |_ 

4- Ab) {9a 9ai) b a 



la 9ai ) i 



(20) 



with g a the actual gauge coupling at scale t > tj (including the extra matter), g ai the gauge 
coupling at the heavy threshold tj, and g a the gauge coupling at scale t > tj in the absence 
of new matter. 

New chiral matter can also modify the running of the soft scalar masses through the S 
parameter. The new matter shifts the running of S by changing the running of g±, but it 
can also contribute to the S term directly at the heavy threshold. Combining these effects, 
the value of S above threshold is 



S(t) 



9i J \9ii 



AS 



(21) 
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S(t) + 



gl-gl 

9i 



Situ) - [^-) A 5. 



where AS is the shift in the value of S at the threshold, gi is the low-scale value of the gauge 
coupling, and S(t) is the value the S term would have in the absence of the new matter. 
Inserting this result into Eq. (JTJ, the effect of the new matter on the running of the soft 
scalar masses through the S term is 



5 A 2 



1 

h 



9i 



(h + Ab) 



9h_ 
9io 



Sn — 



2-\ "I 

9l! 



(h + Ab) 



AS. 



gi 



9U 
9i 



So 



(22) 



As before, the effects of the S term on the soft masses are universal up to the hypercharge 
prefactor. Thus, they still cancel out of the linear combinations given in Eq. (jSJ). 



3.2 Yukawa Effects and Useful Combinations 

We have so far neglected the effects of the MSSM Yukawa couplings on the modified running 
of the soft masses. As a result, the shifts in the running of the soft scalar masses written 
above are family universal. There are also non-universal shifts in the soft scalar masses. 
These arise from the Yukawa-dependent terms in the soft scalar mass beta functions, which 
themselves depend on the Higgs and third-generation soft scalar masses. As a result, the 
low-energy spectrum obtained from a theory with universal scalar masses at the high scale 
and additional GUT multiplets can appear to have non-universal soft masses at the high 
scale if the extra GUT multiplets are not included in the RG evolution. These non-universal 
shifts are usually a sub leading effect, but as we illustrate below they can be significant when 
the super symmetric mass of the new GUT multiplets is within a few orders of magnitude of 
the TeV scale. 

Non-universal mass shifts obscure the relationship between the different MSSM genera- 
tions and the source of supersymmetry breaking. This relationship is closely linked to the 
SUSY flavor problem [30], and possibly also to the origin of the Yukawa hierarchy. For ex- 
ample, third generation soft masses that are significantly different from the first and second 
generation values is one of the predictions of the model of Ref. [31], in which strongly- 
coupled conformal dynamics generates the Yukawa hierarchy and suppresses flavor-mixing 
soft terms. The relative sizes of the high scale soft masses for different families is therefore 
of great theoretical interest. 

Even when there are non-universal shifts from new physics, it is sometimes still possible 
to obtain useful information about the flavor structure of the soft scalar masses at the 
high scale. To a very good approximation, the flavor non-universal contributions to the RG 
evolution of the soft masses are proportional to the third generation Yukawa couplings or the 
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trilinear couplings. There is also good motivation (and it is technically allowed) to keep only 
the trilinear couplings for the third generation. In this approximation, the Yukawa couplings 
and the A terms only appear in the one-loop RG equations for the soft masses through 
three independent linear combinations. Of the seven soft masses whose running depends 
on these combinations, we can therefore extract four mass combinations whose evolution is 
independent of Yukawa effects at one-loop order [32]J^| They are: 



The cancellation in the first two terms occurs because the linear combinations of masses 
correspond to L and B global symmetries. They run only because these would-be symmetries 
are anomalous under SU(2)l and U(l)y The other mass combinations can also be related 
to anomalous global symmetries of the MSSM. 

These mass combinations have the same one-loop RG running as certain combinations of 
masses involving only the first and second generations. For example, the m^ 3 combination 
runs in exactly the same way at one-loop as 



for i = 1,2. If these linear combinations are unequal at the low scale, the corresponding soft 
masses will be non-universal at the high scale. This holds in the MSSM, as well as in the 
presence of any new physics that is flavor universal and respects baryon number. On the 
other hand, m 2 B3 = m 2 B does not imply family-universal high scale masses. For example, 
within a SO (10) GUT, a splitting between the soft masses of the 16's containing the first, 
second, and third generations will not lead to a difference between the mass combinations 
in Eq. (123]) at the low-scale. A similar conclusion holds for the mass combinations m\.. 

In the case of m 2 Xz and my 3 , it is less obvious what to compare them to. The trick here 
is to notice that in the absence of Yukawa couplings, runs in the same way as m 2 Ll since 
they share the same gauge quantum numbers. If the S term vanishes as well, m 2 H also has 
the same RG evolution as m 2 Li . This motivates us to define 



for % = 1,2. These mass combinations can be compared with m 2 Xz and in much the same 
way as for m 2 B . and m\, (although comparing the m^-.'s is only useful for S = 0). They also 
correspond to anomalous global symmetries in the limit that the first and second generation 
Yukawa couplings vanish. 

6 We also assume implicitly that the soft masses are close to diagonal in the super CKM basis, as they 
are quite constrained to be [30) . 




(23) 




(25) 
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The mass combinations listed above can be useful if there is heavy new physics that 
hides the relationships between the high scale masses. For instance, suppose the high scale 
soft spectrum obtained using the MSSM RG equations applied to the measured soft scalar 
masses shows a large splitting between mg 3 and . If the corresponding splitting between 
m B3 and m 2 Bi (at any scale) is very much smaller, this feature suggests that there is new 
physics that should have been included in the RG running, or that there exists a special 
relationship between rrin., rnfj_, and m 2 D . at the high scale. A similar conclusion can be made 
for the other mass combinations. 

3.3 Some Numerical Results 

In our numerical analysis, we follow a similar procedure to the one used in the previous 
section. The MSSM running is performed at two-loop order, and we interface with Sus- 
pect 2.3.4 [10] to compute the low-scale threshold corrections. New physics, in the form of 
vector-like GUT multiplets at an intermediate scale is included only at the one-loop level. 
Unlike the previous section, we define our high-energy spectrum using a simple mSUGRA 
model in the Ab 7^ theory, and include the new physics effects in generating the low-energy 
spectrum. We then evolve this spectrum back up to the unification scale using the MSSM 
RG evolution, with Ab = 0. Our goal is to emulate evolving the MSSM soft parameters in 
the presence of unmeasured and unknown high-scale new physics. 

The running of the soft masses with A^ 5e g = 7 sets of 5 ©5 multiplets with an intermediate 
scale mass of jl = 10 11 GeV is shown in Fig. H]for the mSUGRA parameters itlq = 300 GeV, 
m i/2 = 700 GeV, tan/5 = 10, and A = 0. These parameters are used to find the low energy 
spectrum, which is then RG evolved back up to the high scale with Ab = 0. From this 
figure, we see that a naive MSSM extrapolation of the soft parameters (i.e. with Ab = 0) 
yields predictions for the high-scale soft scalar masses that are significantly larger than the 
actual values. In the same way, the MSSM predicted values of the high-scale gaugino masses 
are smaller than the correct values, as can be seen in Fig. [3j Note that for jl = 10 11 GeV, 
^5ffi5 = 7 is about as large as possible while still keeping the gauge couplings perturbatively 
small all the way up to Mqut- 

Besides confusing the relationship between the high scale gaugino and scalar soft masses, 
heavy GUT multiplets can also obscure the comparison of the high scale scalar masses 
from different generations. As discussed above, this arises from the backreaction in the 
Yukawa- dependent terms in the RG equations for the third generation soft scalar masses. 
Numerically, we find that the splitting is quite small compared to the absolute scale of the 
masses for jl > 10 11 GeV. This is illustrated in Fig. We also find that an approximate 
preservation of universality persists for other values of mo, Aq, and tan f3 as well. The reason 
for this appears to be that for jl > 10 11 GeV, the Yukawa couplings are smaller than the 
gauge couplings by the time the new physics becomes relevant. 

A much greater splitting between the high scale values of and ttiq 3 , and m 2 Jl and 
m^ 3 , is obtained for lower values of jl. This effect is shown in Fig. for A^ 5ffi g = 3 sets 
of 5 © 5 multiplets with an intermediate scale mass of jl = 10 4 GeV, and the mSUGRA 
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Figure 4: Running of the first generation soft scalar masses with N 5 q$ = 7 and jl = 10 11 GeV for 
the mSUGRA input parameters tuq = 300 GeV, m\/2 = 700 GeV, tan/3 = 10, and A$ = 0. The 
dashed lines show the actual running of these parameters, while the solid lines show the running 
from low to high using the RG equations of the MSSM, ignoring the additional heavy multiplets. 



parameters mo = 300 GeV, myj = 700 GeV, tan/3 = 10, and Aq = 0. The value of the 
gauge couplings at unification here is very similar to the p, = 10 11 GeV and AT 5ffi 5 = 7 case. 
As might be expected, the Yukawa-dependent terms in the soft scalar mass RG running 
become important at lower scales where the top Yukawa approaches unity. 

It is also interesting to note that in both Figs. [5] and El the soft masses appear to take 
on family universal values, = and = mfj 3 , at the same scale, near 10 15 GeV in 
Fig. EJ and close to 10 10 GeV in Fig. [Si It is not hard to show, using the mass combinations 
in Eqs. fl23f24f25l) . that this feature holds exactly at one-loop order provided S = 0, the high 
scale masses are family-universal, and the only relevant Yukawa coupling is that of the top 
quark. In this approximation, all the family- dependent mass splittings are proportional to 
(nijj — m^ ), and hence vanish when m 2 H = m 2 L . This relationship can be seen to hold 
approximately in Fig. [HI which also includes two-loop and bottom Yukawa effects. 

In Figs. [7] and [S] we show the running of the mass combinations m Bi and m Xi (where 

mi = m 2 / \/\ui 2 \) for N 5S) ^ = 3 and jl = 10 4 GeV with the high scale mSUGRA input values 
m Q = 300 GeV, M ly / 2 = 700 GeV, tan/3 = 10, and A Q = 0. These figures also show the values 
of m Bi and m Xl that would be obtained by running up without including the effects of 
the heavy new physics. Comparing these figures to Figs. and El it is apparent that the 
splittings between m% and m B3 , and m 2 Xl and rn X3 , are very much less than the high scale 
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log 10 (Q/GeV) 

Figure 5: Running of the soft scalar masses of and U\g with N 59 g = 7 and jl = 10 11 GeV for 
the mSUGRA input parameters tuq = 300 GeV, m\/2 = 700 GeV, tan/3 = 10, and Aq = 0. The 
dashed lines show the actual running of these parameters, while the solid lines show the running 
from low to high using the RG equations of the MSSM, ignoring the additional heavy multiplets. 



splittings between the Q and U soft masses. 

These relationship between the B and X soft mass combinations from different families 
is a footprint left by the full theory (including the heavy GUT multiplets) on the low-energy 
spectrum. Since the scalar masses in the full theory are universal at the high scale, the low 
scale splittings between the B and X soft mass combinations are very small. On the other 
hand, running the low scale Q and U scalar soft masses up within the MSSM does not suggest 
any form of family universality among these masses. Therefore, as we proposed above, the 
low-energy values of these particular combinations of soft scalar masses can provide evidence 
for heavy new physicsEl 

The analysis in this section shows that even though new physics in the form of additional 
heavy GUT multiplets can significantly disrupt the predictions for the high scale soft spec- 
trum obtained by running in the MSSM, certain key properties about the input spectrum can 
still be deduced using collider scale measurements. Most significantly, the low-scale values 
of the gaugino masses and gauge couplings can be used to predict the approximate ratios of 
the high-scale values, provided gauge unification is preserved. 

The effect of additional GUT multiplets on the scalar soft masses is more severe. Extrap- 

7 We have checked that the small splittings between m Bi and w B3 , as well as between m Xi and rn Xs> 
arise from higher loop effects. 
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Figure 6: Running of the soft scalar masses of Qis and £/i 3 with iV 5 ^g = 3 and Jl = 10 GeV for 
the mSUGRA input parameters tuq = 300 GeV, mi/2 = 700 GeV, tan/3 = 10, and Aq = 0. The 
dashed lines show the actual running of these parameters, while the solid lines show the running 
from low to high using the RG equations of the MSSM, ignoring the additional heavy multiplets. 



olating the soft masses without including the contributions from the heavy GUT multiplets 
leads to a prediction for the input soft masses that are generally too low. The splittings 
between the soft masses from different generations can be shifted as well. Despite this, some 
of the flavor properties of the input soft mass spectrum can be deduced by comparing the 
evolution of the mass combinations in Eqs. (|23|) . ([MP , and (1251) . For example, m^ 3 = m? Bi 
and m 2 Li = m| 3 suggests some form of flavor universality (or an embedding in SO (10)), even 
if the scalar masses extrapolated within the MSSM do not converge at M GUT . We expect 
that these special mass combinations could prove useful for studying other types of heavy 
new physics as well. 

4 The (S)Neutrino Connection 

We have seen in the previous sections that taking the unification of the gauge couplings 
as a serious theoretical input still leaves considerable room for experimental uncertainties 
and new physics to modify the extrapolated values of the soft supersymmetry breaking 
parameters at very high energies. For example, a Fayet-Iliopoulos D-term for hypercharge 
or additional complete GUT multiplets with intermediate scale masses will not disrupt gauge 
coupling unification, but will in general change the running of the parameters of the model. 
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Figure 7: Running of the soft scalar mass combinations mg 3 and with iV^g = 3 and fl = 
10 GeV for the mSUGRA input parameters mo = 300 GeV, mjrt = 700 GeV, tan/3 = 10, and 
Aq = 0. The dashed lines show the actual running of these parameters, while the solid lines show 
the running from low to high using the RG equations of the MSSM, ignoring the additional heavy 
multiplets. The small deviations in these figures arise from higher loop effects. 



In this section we wish to study the effect of additional intermediate scale singlet matter 
with significantly large Yukawa couplings to the MSSM matter fields. A particularly well- 
motivated example of this, and the one we consider, are heavy singlet neutrino multiplets. 

The observed neutrino phenomenology can be accommodated by extending the matter 
content of the MSSM to include at least two right handed (RH) neutrino supermultiplets 
that are singlets under the SM gauge group [5] . Throughout the present work, we will assume 
there are three RH neutrino flavors. The superpotential in the lepton sector is then given by 

W t = y e LH d E + y„ LH U N R - ^M n N R N R , (26) 

where L, E, and N R are respectively the SU{2) L doublet, SU{2) L singlet, and neutrino chiral 
supermultiplets, each coming in three families. The quantities y e , y u , and Mr are 3x3 
matrices in lepton family space. The and H u fields represent the usual Higgs multiplets. 
The gauge-invariant interactions among leptons and Higgs superfields are controlled by the 
family-space Yukawa matrices y e and y„. As is conventional, we shall implicitly work in a 
basis where y e is diagonal. Since the N R are singlets, we can also add to the superpotential 
a Majorana mass Mr for these fields. 

Assuming the eigenvalues of Mr lie at a large intermediate mass scale, 10 9 — 10 14 GeV, 
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Figure 8: Running of the soft scalar mass combinations mx 3 and mx 1 with iV^g = 3 and 
p, = 10 4 GeV for the mSUGRA input parameters tuq = 300 GeV, = 700 GeV, tan/3 = 10, and 
Aq = 0. The dashed lines show the actual running of these parameters, while the solid lines show 
the running from low to high using the RG equations of the MSSM, ignoring the additional heavy 
multiplets. The small deviations in these figures arise from higher loop effects. 



we can integrate out the RH neutrino superfields and obtain a term in the effective super- 
potential that leads to small neutrino masses through the see-saw mechanism, 

W m „ = -iyjM^y* LH U LH U . (27) 
After electroweak symmetry breaking, the neutrino mass matrix becomes 

(m i ,)=yjM R 1 y,^, (28) 
where v u = (H u ). The mass matrix can be conveniently diagonalized by the transformation 

(mf*) = U T (m,)U, (29) 

with U a unitary matrix. This matrix U is the usual PMNS matrix that describes lepton 
mixing relative to the flavor basis where the charged lepton Yukawa matrix y e is diagonal. 
One can also write the neutrino Yukawa matrix as [33J 

y^-^M^R ^/m^Ut (30) 



20 



where v u = (H u ) and R is a complex orthogonal matrix that parametrizes our ignorance of 
the neutrino Yukawas. As an estimate, Eq. (1301) shows that the size of the neutrino Yukawa 
couplings will be on the order of 

0-57 f M R \ 3 / m „ \ I 
^ sin/3 Vl0 14 GeVy' WleV/ ' 1 j 

Thus, the neutrino Yukawa couplings can take large 0(1) values comparable to the top 
Yukawa coupling for Mr ~ 10 14 GeV. If this is the case, then above the see-saw mass 
threshold the effects on the RG running of the MSSM soft parameters due to the neutrino 
Yukawas can be substantial. 

The addition of RH neutrinos to the MSSM can lead to lepton flavor violation (LFV) 
through the RG running of the off-diagonal slepton mass terms [3H [35] . In this work we will 
only consider simple scenarios of neutrino phenomenology in which the amount of lepton 
flavor violation induced by the heavy neutrino sector is small. However, the observation of 
LFV signals could potentially provide information about a heavy neutrino sector [3J1 [351 
l36l 137] . Precision measurements of the slepton mass matrices can also be used to constrain 
possible heavy neutrino sectors [38| [39] . Heavy singlet neutrinos may also be related to the 
source of the baryon asymmetry through the mechanism of leptogenesis pTO] . 



4.1 Running Up 

The strategy we use in this section is similar to the one followed in the previous sections. We 
assume a universal high scale mass spectrum at Mqut, and RG evolve the model parameters 
down to the low scale Mi ow = 500 GeV including the additional effects of the neutrino sector 
parameters. The resulting low scale spectrum is then run back up to Mqut using the RG 
equations for the MSSM without including the neutrino sector contributions. As before, we 
use this procedure to illustrate the discrepancy between the extrapolated parameter values 
and their true values if the new physics effects are not included in the running. 

To simplify the analysis, we make a few assumptions about the parameters in the neutrino 
sector. We choose the orthogonal matrix R to be purely real, and we take the heavy neutrino 
mass matrix Mr to be proportional to the unit matrix, Mr = Mr I. We also take the 
physical neutrino masses to be degenerate. This allows us to write 

so that the neutrino Yukawa couplings have the form of a universal constant multiplying a 
unitary matrix. To fix the value of y v , we will set m v = 0.1 eV@. This choice is close to being 

8 The diagonal neutrino mass matrix m^ lag and the lepton mixing matrix U are measured at low scales, 
and one should really evaluate them at the intermediate scale Mr by running the Yukawa couplings up [37] . 
Since we are most interested in the effect of the neutrino Yukawas after reaching the intermediate scale, we 
will neglect this additional running below Mr. 
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as large as possible while remaining consistent with the cosmological bounds on the sum of 
the neutrino masses, < 0.68 eV [H] . Note that larger values of the neutrino masses tend 
to maximize the size of the resulting neutrino Yukawa couplings. 

We also need to impose boundary conditions on the RH sneutrino masses and trilinear 
couplings. We take these new soft parameters to have universal and diagonal boundary 
conditions at the unification scale, 

m2 N l3 = m l5ij, and = A Q y Vij , (33) 

where m and A are the same universal soft scalar mass and trilinear coupling that we 
will apply to the MSSM fields in the analysis to follow. With these choices for the neutrino 
parameters, the effects of the neutrino sector on the (one-loop) RG running of the MSSM 
soft terms take an especially simple form. In particular, the amount of leptonic flavor mixing 
induced is expected to be very small, and the diagonal and universal form of the neutrino 
sector soft terms, Eq. ( 1331) . is approximately maintained at lower scales. 

To illustrate the effects of the neutrino sector, we set the high scale spectrum to coincide 
with the SPS-5 benchmark point, and we extend the corresponding soft terms to the neutrino 
sector. The input values at Mqut for this point are m = 150 GeV, mi/2 = 300 GeV, 
Aq = —1000 GeV, and tan/3 = 5 with /i positive. These input values tend to magnify the 
effects of the neutrino sector because the large value of A Q feeds into the running of the soft 
masses. Large neutrino Yukawa couplings alter the running of the top Yukawa coupling as 
well. 

Of the MSSM soft parameters, the greatest effects of the heavy neutrino sector are seen 
in the soft scalar masses and the trilinear A terms. The gaugino masses are only slightly 
modified. The evolution of the soft masses for the H u , L3, and U3 fields from low to high 
are shown in Fig. [91 both with and without including the effects of the neutrino sector for 
Mr = 3 x 10 14 GeV. For this value of Mr and with tan (3 = 5, the Yukawa coupling is 
close to being as perturbatively large as possible. The extrapolated values of rrijj and m\ 3 
deviate significantly from the actual input values if the effects of the neutrino sector are not 
taken into account in the RG evolution. These fields are particularly affected because they 
couple directly to the heavy neutrino states through the neutrino Yukawa coupling. The 
shift in the running of m^ 3 arises indirectly from the effect of the neutrino Yukawas on m 2 Hu 
and the top Yukawa coupling y t . 

In Fig. [10] we show the size of the discrepancies in the extrapolated high scale values of 
a few of the soft scalar masses if the neutrino sector effects are not included in the running. 
These discrepancies are plotted as a function of the heavy neutrino mass scale Mr. As 
above, the high scale input spectrum consists of the SPS-5 values. Again, the soft masses 
m 2 Hu and m 2 L . are altered the most, although the third generation squark soft masses also 
get shifted somewhat as a backreaction to the changes in m 2 H and the top Yukawa coupling. 
This plot also shows that the sizes of the discrepancies remain quite small for Mr less than 
10 13 GeV. Values of Mr considerably less than this are favored if leptogenesis is to be the 
source of the baryon asymmetry [10]. Neutrino masses well below 0.1 eV would also lead to 
less pronounced deviations in the extrapolated soft masses. 
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Figure 9: Running of the soft scalar masses of H u , L%, and C/3 for SPS-5 input parameters at 
Mqut with three additional heavy RH neutrinos of mass M R = 10 14 GeV. The solid lines show 
the full running, including the neutrino sector effects, while the dashed lines show the low-to-high 
running of the soft masses in the MSSM, with the neutrino sector effects omitted. 



The deviations induced by not including the new neutrino sector physics in the running 
take a similar form to those obtained by not taking account of heavy GUT multiplets. 
For both cases, the gaugino mass running is only modified in a very controlled way, while 
the soft scalar masses and trilinear couplings deviate more unpredictably. In particular, 
the high scale flavor structure of the soft masses can be obscured. With large neutrino 
Yukawa couplings, the third generation squark masses receive additional contributions to 
their running relative to the first and second generations due to the potentially large effect 
of the neutrino Yukawa couplings on m 2 Hu . The sizes of these additional family-dependent 
shifts tend to be fairly small, as can be seen in Fig. [TUJ Furthermore, these effects cancel 
out in the mass combinations m 2 B . denned in Eqs. fl23l) and (l24j) . and we find m Bl ~ m 2 B3 at 
all scales, regardless of whether of not the neutrino effects are included. On the other hand, 
the running of m\ 3 relative to m\ , m ^c 3 relative to m 2 Xv and relative to raf^ need no 
longer coincide if there is a heavy neutrino sector. 

Finally, let us also mention that for more general neutrino sector parameters than those 
we have considered, there can arise significant lepton flavor mixing couplings in the MSSM 
slepton soft terms from the RG running [34]. Measurements of this mixing in lepton fla- 
vor violating processes can therefore provide an experimental probe of the heavy neutrino 
multiplets [31 1 135 1 136 1 137]. A measured splitting among the three slepton masses m\ would 
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Figure 10: The high scale (Mgut) values of the soft scalar masses extrapolated using the MSSM 
RG equations, without including neutrino sector effects. The low-scale (500 GeV) values of the soft 
masses used in the extrapolation were obtained from SPS-5 input parameters at Mgut with three 
additional heavy RH neutrinos of mass Mr. The deviations from rrii(MGUT) = 150 GeV represent 
the discrepancy between the MSSM extrapolated values and the correct value in the full theory 
with heavy RH neutrinos. These discrepancies are shown as a function of the heavy neutrino scale 
Mr. 



also constitute another indication of the existence of a neutrino sector with sizeable Yukawa 
couplings and nontrivial flavor structure. Both high and intermediate energy data may be 
complementary and very useful in extrapolating the MSSM soft terms to high energies. 

5 Putting it All Together: an Example 

In this section we summarize some of our previous results with an explicit example. We 
begin with a low energy spectrum of MSSM soft supersymmetry breaking parameters that 
we assume to have been measured at the LHC to an arbitrarily high precision. We then 
attempt to deduce the essential features of the underlying high scale structure by running 
the low energy parameters up and applying some of the techniques discussed in the previous 
sections. 

In our example we will make the following assumptions: 
• The possible types of new physics beyond the MSSM are: 
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— Complete 5 © 5 GUT multiplets with a common (SUSY) mass scale /L 

— Three families of heavy singlet (RH) neutrinos at the mass scale Mr. 

— A fundamental hypercharge D term. 

In the case of complete GUT multiplets, we will assume further that there are no 
superpotential interactions with the MSSM states as in Section [31 For heavy RH 
neutrinos, we will make the same set of assumptions about the form of the mixing and 
mass matrices as in Section HI 

• The high scale spectrum has the form of a minimal SUGRA model (up to the scalar 
mass shifts due to a hypercharge D-term) at the high scale Mqut — 2.5 x 10 16 GeV. 

• This mSUGRA spectrum also applies to the soft parameters corresponding to any new 
physics sectors. For example, a trilinear A term in the RH neutrino sector has the 
form sl u = Aq y p at Mqut, where Aq is the universal trilinear parameter. 

These assumptions are not entirely realistic, but they make the analysis tractable. More- 
over, even though this exercise is highly simplified compared to what will be necessary should 
the LHC discover supersymmetry, we feel that it illustrates a number of useful techniques 
that could be applied in more general situations. With this set of assumptions, the underlying 
free parameters of the theory are: 

mo, mi/2, Aq, f, A 5e g, fl, Mr. (34) 

where mo, my 2 and Aq are common mSUGRA inputs at Mqut, £ is the fundamental hy- 
percharge D term, A 5e § is the number of additional 5 © 5 multiplets in the theory with a 
super symmetric mass fl, and Mr is the heavy neutrino scale. 

5.1 Step 1: Running Up in the MSSM 

As a first step, we run the low energy spectrum up to the high scale M GUT using the RG 
equations for the MSSM, without including any potential new physics effects. The low 
energy MSSM soft spectrum we consider, defined at the low scale Mi ow = 500 GeV, is 
given in Table HJ In addition to these soft terms, we also assume that tan (3 = 7 has been 
determined, and that the first and second generation soft scalar masses are equal. With 
this set of soft terms, we have verified that the low energy superpartner mass spectrum is 
phenomenologically acceptable using SuSpect 2.3.4 [ID] . The lightest Higgs boson mass is 
114 GeV for a top quark mass of m t = 171.4 GeV. 

Even before extrapolating the soft parameters, it is possible to see a number of interesting 
features in the spectrum. The most obvious is that the low scale gaugino masses have ratios 
close to Mi : M 2 : M 3 ~ 1 : 2 : 6. This suggests that the high scale gaugino masses 
have a universal value my 2, and provides further evidence for gauge unification. The low- 
energy value of the S term, as defined in Eq. §2§, is also non-zero and is in fact quite large, 
S{M low ) ~ (620 GeV) 2 . This indicates that there are significant contributions to the effective 
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hypercharge D term in the high scale theory. Since the S term is non-zero, it is also not 
surprising that Sb-l — (446 GeV) 2 , as defined in Eq. (jTOl) . is non-zero as well. 



Sort Parameter 


T (~1 1 T T I 

Low Scale Value 


TT* 1 d 1 T T 1 

High Scale Value 




(GeV) 


(GeV) 


Mi 


146 


356 


M 2 


274 


355 


M 3 


859 


370 




-956 


-766 




-1755 


-818 




-737 


-524 




-700 


419 




350 


236 


m Qz 


821 


549 


m Ua 


603 


445 


m D3 


884 


501 


m u 


356 


213 


m E3 


349 


404 


m Qi 


934 


532 


m Vx 


872 


402 


m Dl 


888 


501 


m Ll 


357 


213 


m El 


352 


404 



Table 1: The low-energy scale {M\ ow = 500 GeV) soft supersymmetry breaking spectrum used in 
our analysis. The soft scalar masses listed in the table correspond to the signed square roots of the 
actual masses squared. In this table we also the high scale values of these soft parameters obtained 
by running them up to Mqut — 2.5 x 10 16 GeV using the RG evolution appropriate for the MSSM. 

The values of the soft parameters extrapolated to Mqjjt within the MSSM are listed in 
Tabled) The MSSM running of the soft scalar masses is also shown in Fig.[TTJ As anticipated, 
the gaugino masses unify approximately to a value M\ ~ M2 — M3 ~ mi/2 = 350 GeV at 
Mqut- The high scale pattern of the soft scalar masses (and the trilinear A terms) shows 
less structure, and is clearly inconsistent with mSUGRA high scale input values. 

Since S(Mi ow ) is large and non-zero, we are motivated to look for a hypercharge D 
term contribution to the soft scalar masses. Such a contribution would cancel in the mass 
combinations 

Ami = ( Y J m * - Y i m ))/^ ~ Y i)- (35) 

If the high scale soft masses have the form mf = + \J\di Yi £ ; as in mSUGRA with a 
hypercharge D term, these combinations will all be equal to m^ at this scale. The high 
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scale soft masses here, extrapolated within the MSSM, exhibit no such relationship. Even 
so, these mass combinations will prove useful in the analysis to follow. 

It is also interesting to compare the pairs of mass combinations m\., m 2 B ., m 2 ^., and my 
for i = 1, 3, as defined in Section [31 Of these, the most useful pair is m 2 Bl and m^. At the 
low and high scales (extrapolating in the MSSM) this pair obtains the values 

m 2 Bl (M low ) ~ (441 GeV) 2 , m 2 Ba (M low ) ~ (452 GeV) 2 , 
m 2 Bl {M GUT ) ~ (392 GeV) 2 , m 2 B3 {M GUT ) ~ (392 GeV) 2 . {60> 

The near equality of m 2 Bl and m B at the high scale is particularly striking. By comparison, 
there is a significant inter-family splitting that occurs between the high scale soft masses of 
Ui and U 3 , 

ml^Maur) = (402 GeV) 2 , m 2 Uz {M GUT ) = (445 GeV) 2 . (37) 

This apparent fine-tuning among the soft masses that make up the 5-type combinations 
is suggestive of an underlying structure in the theory. However, based on the values of 
the individual soft masses extrapolated within the MSSM, such a structure is not obvious. 
Instead, we can interpret this as a hint for new intermediate scale physics. 

Without our guiding assumptions, the high scale spectrum listed in Table [H obtained by 
running up in the MSSM does not exhibit any particularly remarkable features aside from 
the universality of the gaugino masses. Even so, the curious relationship between the m 2 B 
and m 2 B mass combinations provides a strong hint that we are missing something. It is not 
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clear how strong this hint would have been had we also included reasonable uncertainties in 
the low scale parameter values. 



5.2 Step 2: Adding GUT Multiplets 

As a first attempt to fit the low energy soft spectrum to the class of models outlined above, 
let us consider adding additional vector-like GUT multiplets the the theory at the scale 
fx. We try this first because, as we found in Sections [3] and EJ the contributions from such 
multiplets are potentially much larger than those due to heavy singlet neutrinos. 

In adding the new GUT multiplets, we will make use of our starting assumptions about 
the possible forms of new physics. Given the large value of S(Mi at0 ), there appears to be 
significant hypercharge D term. Also from our assumptions, this D term will contribute to 
the soft scalar masses of the heavy GUT multiplets, which will in term feed into the running 
of the MSSM scalar masses through the S term above the scale fx. For this reason, it is safer 
to work with the mass differences defined in eq. (135]) whose running (to one-loop) does not 
depend on the S-term. 

Among the low scale soft masses listed in Table [fl we expect the slepton soft mass m 2 El 
to be among the easiest to measure, and the least susceptible to new physics effects. Thus, 
we will use it as a reference mass in all but two of the differences we choose. The mass 
differences we consider are 

Am QiBn Am u,Ev Am D lEl > Am L lEl 

Am2 H d EV Am H uEl i Am Q lDl , Am Q 3 D^ ( 38 ) 
Am Q 3 E^ Am U 3 Ev Am D 3 E V Am i 3 £i- 



These depend on several independent mass measurements. 

In Fig. [12] we show the high scale values of these mass differences obtained by running 
up the low scale soft masses while including a given number of iV 5e 5 additional 5 © 5 
GUT multiplets at the scale fx. For each of the plots, nearly all the mass differences unify 
approximately, as they would be expected to do if the underlying theory has a mSUGRA 
spectrum. The best agreement with a mSUGRA model is obtained for N 5 ^ = 5 with 
fx = 10 10 GeV. (More precisely, the agreement is obtained when a shift Ab = —5 is applied 
to the gauge beta function coefficients 6j at the scale fx = 10 10 GeV). 

Taking N 5(B g = 5 and fx = 10 10 GeV, the high scale values of the mass differences are 

Am 2 QiEl = (200 GeV) 2 = Am 2 UiEi = Am 2 DlEi = Am 2 HdEl (39) 

Am 2 QaEi = (197 GeV) 2 , Am 2 ^ = (200 GeV) 2 
Am^ £l = -(157 GeV) 2 , Am 2 ^ = (183 GeV) 2 

Most of these high scale values coincide, suggesting a mSUGRA value for the universal soft 
scalar mass of about mo = 200 GeV. On the other hand, the soft mass differences involving 
the H u and L fields show a significant deviation from this near-universal value. Based on 
the results of Section HJ these are precisely the scalar masses that are most sensitive to a 
heavy singlet neutrino sector. 
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Using this same choice of new physics parameters, we can also estimate the values of the 
other mSUGRA parameters. Heavy singlet neutrinos are not expected to significantly alter 
the running of the gaugino soft mass parameters. If we run these up to Mqut including 
-W505 — 5 additional GUT multiplets at jl — 10 10 GeV, we find mi/ 2 — 700 GeV. Doing the 
same for the trilinear couplings, we do not find a unified high scale value for them. Instead, 
we obtain A t (M GUT ) = -401 GeV, A T = -407 GeV, and A b = -500 GeV. This is not 
surprising since a heavy RH neutrino sector would be expected to primarily modify A t and 
A T , while having very little effect on A^. Thus, we also expect A ~ —500 GeV. 



> 

a 




H E. 

u 1 




7 8 9 9 10 11 10 11 12 

log . n (£/GeV) log. n (fL/GeV) log . .(ft /GeV) 



-le+05 



Figure 12: 
function of 



High scale values for the soft scalar mass differences Am| defined in Eq. (|38p as a 
the mass scale jl of the iV^g heavy GUT multiplets for = 3, 5, 7. 



It is possible to estimate the value of the hypercharge D term as well. Using the hypoth- 
esis mf = mQ + \J\gi Y£ at the high scale, we find 

£(AW) = J\- (m 2 El - m 2 Hd )/(Y E - Y Hd ) ~ (494 GeV) 2 . (40) 

We obtain similar values from the corresponding combinations of other mass pairs with the 
exception of L and H u . Based on our previous findings, we suspect that the L and H u soft 
masses are modified by a heavy RH neutrino sector. 
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Note that had we included experimental and theoretical uncertainties it would have been 
considerably more difficult to distinguish different values of A^g and jl. Instead of finding 
a single value for 5 © 5 and a precise value for jl, it is likely that we would have only been 
able to confine A^g and \x to within finite ranges. 

5.3 Step 3: Adding a Heavy Neutrino Sector 

By adding iV 50 g = 5 complete 5 © 5 multiplets at the scale jl = 10 10 GeV and a hypercharge 
D term, we are nearly able to fit the low scale spectrum given in Table [1] to a mSUGRA 
model with mo = 200 GeV, = 700 GeV, and Aq = —500 GeV. However, there are 

several small deviations from this picture, most notably in the soft masses for H u and L as 
well as the trilinear couplings A t and A T . We attempt to fix these remaining discrepancies 
by including heavy RH neutrinos at the scale Mr. 

Given our initial assumptions about the form of a possible RH neutrino sector, the only 
independent parameter in this sector is the heavy mass scale Mr. To investigate the effects of 
RH neutrinos, we examine the high scale values of the mass differences given in Eq. (135]) for 
the third generation scalars and H u , using m El as a reference mass. We add a RH neutrino 
sector with heavy mass Mr and run the low scale parameters listed in Table [I] subject to 
the additional neutrino effects, as well as those from Af 5e g = 5 heavy GUT multiplets with 
jl = 10 1 GeV. The result is shown in Fig. [13], in which we scan over Mr. This plot shows 
that if we include heavy RH neutrinos at the mass scale near Mr = 10 14 GeV, all the mass 
differences will flow to a universal value of about mo = 200 GeV at the high scale. 

We can further confirm this result by examining the high scale trilinear couplings obtained 
by this procedure. These also attain a universal value, A = —500 GeV, at the high scale 
for Mr = 10 14 GeV as shown in Figure [TH These universal values are consistent with those 
we hypothesized before the inclusion of heavy neutrino sector effects. Similarly, we can also 
study the value of £ obtained using Eq. fHOj) . but using the high scale L and H u soft masses 
computed by including heavy RH neutrinos in their RG evolution. As before, we obtain a 
value f ~ (494 GeV) 2 . 

5.4 Summary 

We have succeeded in deducing a high scale mSUGRA model augmented by heavy new 

physics effects that reproduces the soft spectrum in Table [TJ The relevant parameter values 
are: 

m = 200 GeV, m 1/2 = 700 GeV, A = -500 GeV, 

AW = 5, A = 10 10 GeV, (41) 

Mr = 10 14 GeV, £ = (494 GeV) 2 . 

Our example did not include potential uncertainties in the input soft parameter values. It 
is likely that such uncertainties would make the analysis more challenging. 

Running the low scale soft parameters up to M GUT within the MSSM, and without 
including any new physics, we obtained a reasonable but mostly unremarkable high scale 
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Figure 13: High scale values for third family and Higgs boson soft scalar mass differences Am| 
(as defined in Eq. (|35p ) as a function of the mass scale Mr of the heavy RH neutrinos. ^5^5 = 5 
GUT multiplets were included in the running above the scale fx = 10 10 GeV. 



soft spectrum. The most obvious feature of this spectrum is the unification of the gaugino 
masses. A more subtle aspect of the high scale spectrum is the small splitting between 
m 2 Bl and m 2 B , A relative to that between the and rriQ 3 , and and mfj 3 . This feature 
hinted at an underlying family-universal flavor spectrum obscured by new physics effects. 
It is not clear whether this hint would survive in a more complete treatment that included 
uncertainties in the input parameter values. By adding new physics, in the form of heavy 
GUT multiplets and RH neutrinos, a simple mSUGRA structure emerged. 

6 Conclusions 

If supersymmetry is discovered at the LHC, the primary challenge in theoretical particle 
physics will be to deduce the source of supersymmetry breaking. By doing so, we may 
perhaps learn about the more fundamental theory underlying this source. In most models of 
supersymmetry breaking, the relevant dynamics take place at energies much larger than those 
that will be directly probed by the LHC. It is therefore likely that the soft supersymmetry 
breaking parameters measured by experiment will have to be extrapolated to higher scales 
using the renormalization group. Given the apparent unification of gauge couplings in the 
MSSM only slightly below Mpi, we may hope that there is little to no new physics between 
the LHC scale and the supersymmetry breaking scale so that such an extrapolation can be 
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Figure 14: High scale values for the soft trilinear A terms as a function of the mass scale of the 
heavy RH neutrinos. = 5 GUT multiplets were included in the running above the scale 

fi = 10 10 GeV. 



performed in a straightforward way. 

Gauge unification still allows for some types of new physics at intermediate scales such 
as complete GUT multiplets and gauge singlets. If this new physics is present, RG evolving 
the MSSM soft parameters without including the new physics effects can lead to an incorrect 
spectrum of soft parameters at the high scale. Even without new intermediate physics, if 
some of the MSSM soft parameters are only poorly determined at the LHC, or not measured 
at all, there can arise significant uncertainties in the RG running of the soft masses that 
have been discovered. 

In the present work we have investigated both of these potential obstacles to running up 
in the MSSM. The soft scalar masses are particularly sensitive to these effects, but we find 
that the gaugino soft masses, and their ratios in particular, are considerably more robust. If 
any one of the scalar soft masses goes unmeasured at the LHC, the running of the remaining 
these soft terms can be significantly modified by the effects of the hypercharge S term. These 
effects are especially severe for the slepton soft masses, which otherwise do not tend to run 
very strongly at all. We find that the uncertainties due to the S term can be avoided if we 
consider the soft mass differences (YjTn? — i^m|), where Yi denotes the hypercharge of the 
corresponding field. If all the soft mass are measured, the soft scalar mass combinations S 
and Sb-l, denned in Eqs. fl2]) and ( fTUl) . provide useful information about potential GUT 
embeddings of the theory. 
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We have also investigated the effects of two plausible types of new physics beyond the 
MSSM that preserve consistency with gauge unification; namely complete vector-like GUT 
multiplets and heavy singlet neutrinos. In the case of complete GUT multiplets, extrapolat- 
ing the measured low-energy soft parameter values without including the additional charged 
matter in the RG running leads to high scale predictions for the gaugino masses that are too 
small, and soft scalar masses that are too large. Even so, the ratios of the gaugino masses at 
the high scale are not modified at leading order, and can be predicted from the low-energy 
measured values provided gauge unification occurs. The extrapolated values of the scalar 
masses are shifted in more complicated ways, and relationships such as family universality at 
the high-scale can be obscured. Despite this, certain hints about the underlying flavor struc- 
ture of the soft masses can still be deduced from the properties of special linear combinations 
of the soft masses, such as {2m 2 n 3 — mfj 3 — m 2 D3 ) relative to (2777,^ — — m 2 Di ). 

The running of the MSSM soft masses can also be modified if there are heavy singlet 
neutrino chiral multiplets in the theory. These can induce small masses for the standard 
model neutrinos through the see-saw mechanism. If the singlet neutrino scale is very heavy, 
greater than about 10 13 GeV, the corresponding neutrino Yukawa couplings can be large 
enough to have a significant effect on the running of the soft masses of H u and L. We have 
studied the size of these effects, as well as the shifts in the other soft masses. The extrapolated 
values of the gaugino masses and the squark soft masses are only weakly modified by heavy 
neutrino sector effects. 

In Section [5] we applied the methods described above to a specific example. In this 
example, the scalar masses have a common high scale value, up to a hypercharge D term. 
However, because of the presence of heavy new physics, this simple structure does not emerge 
when the low scale soft scalar masses are extrapolated up to Mqut using the RG equations 
of the MSSM. Based on the low-energy spectrum alone, we were able to deduce the presence 
of the hypercharge D term. The presence of additional new physics was suggested by the fact 
that the splitting between and was considerably larger than the related splitting 
between m 2 B3 and m 2 Bi . By studying the scalar mass combinations of Eq. fl35l) and including 
heavy GUT multiplets and a right-handed neutrino sector, we were able to reproduce the 
low-energy soft spectrum with an underlying mSUGRA model. 

One can also invert this perspective of overcoming new physics obstacles, and instead 
view these obstacles as providing information and opportunities. As the example of Section [5] 
illustrates, analyses of the kind we consider here probe new physics in indirect ways that can 
lead to convincing arguments for its existence or absence. Such analyses may be the main 
way to learn about new physics that cannot be studied directly. 

If the LHC discovers new physics beyond the standard model, it will be a challenge 
to extract the Lagrangian parameters from the data. It may also be difficult to correctly 
extrapolate these parameters to higher scales in order to deduce the underlying theory that 
gives rise to the low energy Lagrangian. In the present work we have begun to study this 
second aspect of the so-called LHC inverse problem, and we have found a few techniques to 
address some of the potential obstacles. However, our study is only a beginning. We expect 
that a number of additional techniques for running up could be discovered with more work. 
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A similar set of techniques could also be applied to understanding the high scale origin of 
other types of new physics beyond the standard model. These techniques deserve further 
study 
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Appendix: Useful Combinations of Scalar Masses 

In this appendix we collect some combinations of soft scalar masses that are particularly 
useful for running up. 

S Term Effects 

The mass differences 

Ami = ( Y i m i ~ Y i m j)/( Y i ~ Y i) ( 42 ) 
are useful when there is a non-vanishing hypercharge D term. A hypercharge D term can 
shift the low energy values of the soft scalar masses, and can also modify their RG running 
through the S term, as discussed in Section [21 and in Refs. [251 [261 H2]- The effects of a 
hypercharge D term cancel out these mass differences, as well as in the RG equations for 
them. 

This feature is helpful for running up because the effect of the D term on the running 
is determined by the low scale value of the S'-term, which depends on all the soft scalar 
masses (of hypercharged fields) in the theory. If one of these soft masses goes unmeasured, 
there will be a large uncertainty in the value of the S term, and this in turn will induce a 
significant uncertainty in the high scale values of the soft scalar masses after RG evolution. 
By focusing on the mass differences of Eq. (142j) . these ambiguities cancel each other out. 

On the other hand, if all the MSSM soft scalar masses are measured, the low scale values 
of the soft scalar mass combinations 

S = Tr(Ym 2 ), (43) 
S B - L = Tr[(B-L)m% 

provide useful information about the high scale theory, and can be used to test possible GUT 
embeddings of the MSSM [26]. 
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Flavor Splitting Effects 

New physics can obscure the underlying flavor structure of the soft scalar masses. Family- 
universal soft masses derived from a theory containing new physics can generate a low-energy 
spectrum that does not appear to be family-universal after it is evolved back up to the high 
scale without including this new physics. This is true even if the new physics couples in a 
flavor universal way to the MSSM. We presented a particular example of this in Section [31 
where the new physics took the form of complete GUT multiplets having no superpotential 
couplings with the MSSM sector. 

There are four pairs of soft mass combinations that are helpful in this regard [32J. By 
comparing these pairs (at any given scale), it is sometimes possible to obtain clues about 
the underlying flavor structure of the MSSM soft masses. These combinations are: 



m % = 2m l 3 - m h 

m B 3 = 2m Qs- m U 3 -™D 3 - m B, = 2^1^-771^-771^ 

m x 3 = 2m H u -3 m u 3 

m\ z = 3 m|, 3 + 2 77l| 3 





m\ 


= 2m 2 Ll 




™ 2 d 3 «-» 


m 2 Bl 


= 2774 


- 777^ - 


<— > 


m\ x 


= 2m| x 




— 2 m 2 Hd <-» 




= 3777^ 





(44) 



If the high scale soft scalar masses are family-universal, we expect that each of these pairs, 
with the possible exception of the m\., to be roughly equal at the low scale in the MSSM. 
The m 2 x . combinations are expected to match only if 5* = as well. 

To apply the soft mass combinations in Eq. (jH]), one should compare them to the splitting 
between individual soft masses after running all soft masses up to the high scale within the 
MSSM (without new physics). For instance, an inequality of the form 

Wb 3 -TObJ < max{\m 2 Q:i -mjj, \m 2 Uz -m^J, \m 2 D:j - m 2 Dl \] , (45) 

is suggestive of high scale family-universality or a particular relationship between the Q, 
U, and D soft masses that has been obscured by new physics. This can arise from GUT 
multiplets as in Section [31 or from a heavy RH neutrino sector as in Section HI Note that 
heavy neutrinos can disrupt the relationships between the A, X, and Y pairs. 



35 



References 



[1] For reviews, see: H. E. Ha ber and G. L. K ane, Phys. Rept. 117, 75 (1985), S. P. Martin, 
A supersymmetry primer, [hep-p h/9709356; D. J. H. Chung, L. L. Everett, G. L. Kane, 
S. F. King, J. D. Lykken and L . T. Wang, Phys. Rept. 407, 1 (2005) |hep-ph/0312378| ; 
M. A. Luty, [hep-th/0509029| . 

[2] S. Dimopoulos and H. Georgi, Nucl. Phys. B 193, 150 (1981); S. Dimopoulos, S. Raby 
and F. Wilczek, Phys. Rev. D 24, 1681 (1981). 

[3] For recent work on this subject, see: P. M. Zerwas et al, |hep-ph/0211076[ B. C. Al- 
lana ch, G. A. Blair, S . Kraml, H. U. Martyn, G. Polesello, W. Porod and P. M. Zer- 
was, hep-ph/0403133l G. Weiglein et al. [LHC/LC Study Group], Phys. Rept. 426, 



47 (2006) |h ep-ph/0410364 ]; R. Lafaye, T. Plehn and D. Zerwas, |hep-ph/0404"282j 
P. Bechtle, K. Desch and P. Wienemann, Comput. Phys. Commun. 174, 47 (2006) 
|hep-ph / 04 1 2~0T2] ; N . Arkani-Hame d, G. L. Kane, J. Thaler and L. T. Wang, JHEP 
06 08, 070 (2006) |hep-ph/05 12190]. 



[4] Y. Yamada, Phys. Rev. D 54, 1150 (1996) |hep-ph/9602279|; D. M. Pierce, J. A. Bag- 
ger, K. T. Matchev and R. j. Zhang, Nucl Phys. B 491, 3 (1997) |hep-ph/9606211| . 
Y. Yamada, Phys. Lett. B 623, 104 (20 05) |hep-ph/0506262] ; S. P. Martin, Phys. Rev. 
D 74, 075009 (2006) |hep-ph/0608026] . 

[5] G . A. Blair, W . Porod and P. M. Zerwas, Phys. Rev. D 63, 017703 (2001) 
|hep-ph/0007T07]. G. A Blair, W. Porod and P. M. Zerwas, Eur. Phys. J. C 27, 263 
(2003) |hep- ph/0210058 ]. G. A. Blair, A. Freitas, H. U. M artyn, G. Pol esello, W. Porod 
and P. M. Zerwas, Acta Phys. Polon. B 36, 3445 (2005) [ hep-ph/0512"084] ; 

[6] S. P. Martin, in Proc. of the APS/DPF/DPB Summer Study on the Future of Particle 
Physics (Snowmass 2001) ed. N. Graf, In the Proceedings of APS / DPF / DPB Sum- 
mer Study on the Future of Particle Physics (Snowmass 2001), Snowmass, Colorado, 
30 Jun - 21 Jul 2001, pp P327. 

[7] M. Carena, P. H. Chankowski, M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. 
Phys. B 491, 103 (1997) lhep-ph/96 12261] . 



[8] For a review, see: R. N. Mohapatra et al., [hep-ph/0510213 



[9] H. Baer, M. A . Diaz, P. Quintana and X. Tata, JHEP 0004, 016 (2000) 
|hep-ph/0002245] . 



[10] A. Djouadi, J. L. Kneur and G. Moultaka, |hep-p h/0211331|. 

[11] M. Dine and A. E. Nelson, Phys. Rev. D 48, 1277 (1993) |hep-ph/9303230|; M. Dine, 
A. E. Nelson and Y. Shirman, Phys. Rev. D 51, 1362 (1995) |he p-ph/9408"384|; M. Dine , 
A. E. Nelson, Y. Nir and Y. Shirman, Phys. Rev. D 53, 2658 (1996) | hep-ph/9507378| . 



36 



[12] For a review, se e G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 (1999) 
[hep-ph/980127l] . 



hep-th/0105137 



hep-th/0111231 



[13] M. A. Luty and R. Sundrum, Phys. Rev. D 65, 066004 (2002) 
M. Luty and R. Sundrum, Phys. Rev. D 67, 045007 (2003) 
R. Harnik, H. Murayama and A. Pierce, JHEP 020 8, 034 ( 2002) | hep-ph/0204l22] 
R. Sundrum, Phys. Rev. D 71, 085003 (2005) |hep-th/04060T2l ; M. Ibe, K. I. Izawa, 
Y. Nakayama, Y . Shinbara and T. Yanagida, Phys. Rev. D 73, 015004 (2006) 
|hep-ph/0506023|; M. Ibe, K. I. Izawa, Y. Nakayama, Y. Shinbara and T. Yanagida, 
Phys. Rev. D 73, 03501 2 (2006) [hep-ph/0509229| ; M. Schmaltz and R. Sundrum, 
JHEP 0611, Oil (2006) [hep-th/060805lj ~ 

[14] M. Dine, P. J. Fox, E. Gorbatov, Y. Shadmi, Y. Shirman and S. D. Thomas, Phys. 
Rev. D 70, 045023 (2004) |hep-ph/0405159] . 



[15] A. G. Cohen, T. S. Roy and M. Schmaltz, |hep-ph/0612100l 



[16] S. P. Martin and M. T. Vaughn, Phys. Rev. D 50, 2282 (1994) |hep-ph/9311~340~| . 

[17] I. Jack and D. R. T. Jones, Phys. Lett. B 473, 102 (2000) |hep-ph/9911491]; I. Jack , 
D. R. T. Jones and S. Parsons, Phys. Rev. D 62, 125022 (2000) |h ep-ph/0007291| ; 
I. Jack and D. R. T. Jones, Phys. Rev. D 63, 075010 (2001) |hep-ph/00103Qll . 

[18] B. C. Allanach et al, in Proc. of the APS/DPF/DPB Summer Study on the Future of 
Particle Physics (Snowm ass 2001) e d. N. Graf, Eur. Phys. J. C 25, 113 (2002) [eConf 
C010630, P125 (2001)] |hep-ph/0202233| . 

[19] H. Ba er, C. h Chen, F. Paige and X. Tata, Phys. Rev. D 53, 6241 (1996) 
[hep-ph/9512383] ; H. Ba er, C. h. Ch en, M. Drees, F. Paige and X. Tata, Phys. Rev. D 
59, 055014 (1999) [hep-ph/9809~223]; H. Ba chacou, I. Hinchliffe and F. E. Paige, Phys. 
Rev. D 62, 015009 (2000) |hep-ph/9907518|; I. Hinchliffe and F. E. Paige, Phys. Rev. 
D 61, 095011 (2000) |hep-ph/9907519] . 

[20] D. Denegri, W. Majerotto and L. Rurua, Phys. Rev. D 60, 035008 (1999) 
|hep-ph/990123T] . 

[21] K. A. Assamagan, J. Guasch, S. Moretti and S. Penaranda, Czech. J. Phys. 55, 
B787 (2005) l hep-ph/0409~l89~| ; R. Kinnunen, S. Lehti, F. Moo rtgat, A. Nikitenko 
and M. Spira, Eur. Phys. J. C 40N5, 23 (2005) lhep-ph/0503075] . 

[22] For reviews, see: J. F. Gunion, L. Poggioli, R. Van Kooten, C. Kao and P. Rowson, 
eConf C960625, LTH092 (1996) |hep-ph/9703330~] ; K. A. Assam agan, Y. Coadou and 
A. Deandrea, Eur. Phys. J. directC 4, 9 (2002) [hep-ph/020312l]; M. Carena and 
H. E. Haber, Prog. Part. Nucl. Phys. 50, 63 (2003 ) |hep-ph/02 08209|; V. Buscher and 
K. Jakobs, Int. J. Mod. Phys. A 20, 2523 (2005) |hep-ph/0504099|. 



37 



[23] A. de Gouvea, A . Friedland and H. Murayama, Phys. Rev. D 59, 095008 (1999) 
[hep-ph/980348l] . 

[24] T. Falk, Phys. Lett. B 456, 171 (1999) |hep-ph/9902352] . 

[25] M. Drees, Phys. Lett. B 181, 279 (1986); J. S. Hagelin and S. Kelley, Nucl. Phys. 
B 342, 95 (1990); Y. Kawamura, H. Murayama and M. Yamaguchi, Phys. Lett. B 
324, 52 (1994) | hep-ph/9402254|. Y. Kawa mura, H. Murayama and M. Yamaguchi, 
Phys. Rev. D 51, 1337 ( 1995) [hep-ph/9406245| ; H. C. Cheng and L. J. Hall, Phys. 
Rev. D 51, 5289 (1 995) |hep-ph/9411276] ; C. F~ Kolda and S. P. Martin, Phys. Rev. 
D 53, 3871 (1996) |hep -ph/9503445l ; S. P. Martin, Phys. Rev. D 61, 035004 (2000) 
|hep-ph/9907550] . 

[26] A. E. Faraggi, J. S. Hagelin, S. Kelley and D. V. Nanopoulos, Phys Rev. D 45, 3272 
(1992); A. Lleyda and C. Munoz, Phys. Lett. B 317, 82 (1993) |hep-ph/9308208| ; 
Y. Kawamura, T. Kobayashi and H. Shimabukuro, Phys. Lett. B 436, 108 (1998) 
|hep-ph/98053~3~6~] ; H. Baer, M. A. Diaz, J. Ferrandis and X. Tata, Phys. Rev. D 
61, 111701 (2000 ) [hep-ph/99072lT] ; M. R. Ramage, Nucl. Phys. B 720, 137 (2005) 
|hep-ph/0412153] . 

[27] A. H. Chamsed dine and H. K. Dreiner, Nucl. Phys. B 447, 195 (1995) 
|hep-ph/9503454] . 

[28] S. P. Martin and P. Ramond, Phys. Rev. D 51, 6515 (1995) |hep-ph/9501244| . 

[29] G. L. Kane, J. D Lykken, B. D. Nelson and L. T. Wang, Phys. Lett. B 551, 146 (2003) 
|hep-ph/0207168] . 

[30] See, for example, F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. 
B 477, 321 (1996) |hep-ph/9604"3~87| j; M. Misiak , S. Pokors ki and J. Rosiek, Adv. Ser. 
Direct. High Energy Phys. 15, 795 (1998) |hep-ph/9703442] . 

[31] A. E. Nelson and M. J. Strassler, JHEP 0009, 030 (2000) |hep-ph/0006~25ll ; 
T. Kobayashi and H. Terao, Phys. Rev. D 64, 075003 (2001) [hep-ph/0103028l ; 
A. E. Nelson and M. J. Strassler, JHEP 0207, 021 (2002) |hep-ph/0104051| . 

[32] L. E. Ibanez and C. Lopez, Nucl. Phys. B 233, 511 (1984); L. E. Ibanez, C. Lopez 
and C. Munoz, Nucl. Phys. B 256, 218 (1985); M. Care na, M. Olechowsk i, S. Pokorski 
and C. E. M. Wagner, Nucl. Phys. B 419, 213 (1994) [hep-ph/93li222] . M. Carena, 
M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. B 426, 269 (1994) 
|hep-ph/9402253] . 

[33] J. A. Casas and A. Ibarra, Nucl. Phys. B 618, 171 (2001) |hep-ph/0103065] . 



38 



[34] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57, 961 (1986). J. Hisano, T. Mo- 
roi , K. Tobe, M . Yamaguchi and T. Yanagida, Phys. Lett. B 357, 579 (1995) 
|hep-ph/ 9501407| ; J. Hisano, T. Moroi, K. Tobe and M. Yamaguchi, Phys. Rev. D 
53, 2442 (1996) [hep-ph/9510309] . 

[35] J. A. Casas, J. R. Espinosa, A. Ibarra and I. Navarro, Nucl. Phys. B 573, 652 (2000) 
|hep-ph/9910420] . 

[36] F. Deppisch, H. Pas, A. Redelbach, R. Ruckl and Y. Shimizu, Eur. Phys. J. C 28, 365 
(2003) |hep-ph/0206122] ; T. Blazek and S. F. King, Nucl. Phys. B 662, 359 (2003) 
|hep-ph/0211368|; A. Masiero, S. K. Vempati and O. Vives, Nucl. Phys. Proc. Suppl. 
137, 156 (200 4) [hep-ph/0405017] ; E. Arganda and M. J. Herrero, Phys. Rev. D 73, 
055003 (2006) [ hep-ph/ 05 104051 

[37] S. T. Petcov, T . Shindou and Y. Takanishi, Nucl. Phys. B 738, 219 (2006) 
[hep-ph/0508243] . 

[38] H. Baer, C. Bala zs, J. K. Mizukoshi and X. Tata, Phys. Rev. D 63, 055011 (2001) 
|hep-ph/0010068] . 

[39] S. Davidson and A. Ibarra, JHEP 0109, 013 (2001) | hep-ph/ 01 04076] . 

[40] A recent review of this topic can be found in: W. Buchmuller, R. D. Peccei and 
T. Yanagida, Ann. Rev. Nucl. Part. Sci. 55, 311 (2005) [hep-ph/ 05 02169] . 

[41] D. N. Spergel et al, |astro-ph/060"3449] . 

[42] This effect can play an important role in models of gaugino mediation of supersymme- 
try breaking. See, for example, 

D. E. Kaplan, G. D. Kribs and M. Schmaltz, Phys. Rev. D 62, 035010 (2000) 
|hep-ph/991l"293]; Z. Chack o, M. A. Luty, A. E. Nelson and E. Ponton, JHEP 0001, 
003 (200 0) [hep-ph/9911323] ; D. E. Kaplan and T. M. P. Tait, JHEP 0006, 020 
(2000) |hep-ph/0 004200|; M. Schmaltz and W. Skiba, Phys. Rev. D 62, 095005 
(2000) jhep-ph/0001172 ]; M. Schmaltz and W. Skiba, Phys. Rev. D 62, 095004 
(2000) [hep^p h/0004210]; H. Baer, C. Balazs, A. Belyaev , R. Dermisek, A. Mafi and 
A. Mustafayev, JHEP 0205 , 061 (2002) [hep-ph/0204108] ; C. Balazs and R. Dermisek, 
JHEP 0306, 024 (2003) || hep-ph/0303161]; W. Buch muller, J. Kersten and K. Schmidt- 
Hoberg, JHEP 0602, 069 (2006) ]hep-ph/0512152] ; W. Buch muller, L Covi, J. Ker- 
sten and K. Schmidt-Hoberg, JC AP 0611, 007 (20 06) [hep-ph/0609142] ; J. L. Evans, 
D. E. Morrissey and J. D. Wells, | hep-ph/ 061 1185] . 



39 



